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Cognitive  Code-Division  Channelization 

Kanke  Gao.  Stella  N.  Batalama,  Member,  IEEE,  Dimitris  A.  Pados,  Member,  IEEE, 
and  John  D.  Matyjas,  Member,  IEEE 


Abstract — We  consider  the  problem  of  simultaneous  power  and 
code-channel  allocation  for  a  secondary  transmitter/receiver  pair 
coexisting  with  a  primary  code-division  multiple-access  (CDMA) 
system.  Our  objective  is  to  find  the  optimum  transmitting  power 
and  code  sequence  of  the  secondary  channel  that  maximize  the 
signal-to-interference-plus-noise  ratio  (SINR)  at  the  output  of 
the  maximum  SINR  linear  receiver,  while  at  the  same  time 
the  SINR  of  all  primary  channels  at  the  output  of  their  max- 
SINR  receiver  is  maintained  above  a  certain  threshold.  This  is  a 
non-convex  NP-hard  optimization  problem.  We  propose  a  novel 
feasible  suboptimum  solution  using  semidefinite  programming. 
Simulation  studies  illustrate  the  theoretical  developments. 

Index  Terms — Code-channel  allocation,  code-division  multiple- 
access,  cognitive  radio,  power  allocation,  semidefinite  program¬ 
ming,  signal-to-interference-plus-noise  ratio. 


I.  Introduction 

RECENT  experimental  studies  [  1]  demonstrated  that  much 
of  the  licensed  radio  spectrum  experiences  low  utiliza¬ 
tion.  Cognitive  radio  (CR)  [2]  emerges  as  a  promising  tech¬ 
nology  to  improve  spectrum  utilization  by  allowing  secondary 
users/networks  to  share  spectrum  licensed  by  primary  users. 
As  licensees,  the  primary  users  are  to  have  guaranteed  access 
to  the  spectrum  [3].  Therefore,  the  underlying  challenge  of  CR 
technology  is  to  ensure  the  Quality-of-Service  (QoS)  require¬ 
ments  of  the  primary  users  and,  simultaneously,  maximize  in 
a  best-effort  context  the  QoS  of  the  secondary  users  [4]-[7]. 

Herein,  we  are  particularly  interested  in  cognitive  radio 
built  around  a  code-division  multiple-access  (CDMA)  primary 
system.  In  contrast  to  frequency  or  time  division  operation 
where  cognitive  secondary  users  may  transmit  opportunis¬ 
tically  in  sensed  spectrum  holes/void  only,  cognitive  code¬ 
division  users  may  in  principle  operate  in  parallel  in  frequency 
and  time  to  a  primary  system  as  long  as  the  induced  spread- 
spectrum  interference  remains  below  a  pre-defined  acceptable 
threshold1.  Power  control  for  cognitive  code-division  systems 
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1  While  early  standardization  and  regulation  discussions  have  begun  [8],  no 
conclusive  “interference  temperature”  rules  and  agreements  have  been  reached 
yet. 


was  considered  in  [9]  under  an  “interference  temperature"  con¬ 
straint  (total  secondary  user  disturbance  power  over  primary 
band).  No  optimization  was  carried  out  with  respect  to  the 
code  channels  (signatures)  of  the  secondary  users  in  [9].  In 
contrast,  in  [10]  a  secondary  code  assignment  scheme  was 
presented  to  minimize  the  mean-square  crosscorrelation  of  the 
secondary  code  with  the  primary  received  signal.  Extension  to 
multiple  secondary  users  was  also  considered  in  the  form  of 
iterative  secondary  code  set  construction.  Under  interference 
minimizing  code  assignments,  bit  rate  and  spreading  factor 
adjustments  for  a  secondary  CDMA  system  were  considered 
in  [11].  Interesting  work  outside  the  framework  of  CDMA  CR 
in  the  form  of  joint  beamforming  and  power  allocation  algo¬ 
rithms  was  reported  in  [12],  [13],  while  auction  mechanisms 
for  power  control  were  presented  in  [14]. 

In  this  paper,  we  consider  the  opportunity  to  establish  a 
secondary  code-division  link  coexisting  with  a  primary  CDMA 
system.  In  particular,  we  study  for  the  first  time  the  problem  of 
designing  a  joint  power  and  code-channel  allocation  scheme 
for  the  secondary  link  that  maximizes  the  output  Signal-to- 
Interference-plus-Noise-Ratio  (SINR)  of  the  maximum-SINR 
linear  receiver  filter  under  SINR  constraints  on  all  primary 
users2  and  a  peak  transmission  power  constraint  on  the 
secondary  user.  We  recognize  that,  regretfully,  this  key  CR 
CDMA  formulation  is  a  non-convex  NP-hard  problem.  Yet, 
using  semidefinite  programming  methodology  we  are  able  to 
develop  a  novel,  realizable  suboptimum  solution  with  excellent 
cognitive  system  performance  characteristics  as  demonstrated 
by  simulation  studies  included  in  this  paper. 

The  rest  of  the  paper  is  organized  as  follows.  Section 
II  is  devoted  to  the  description  of  the  CR  CDMA  system 
model  and  our  formulation  of  the  optimization  problem.  In 
Section  III,  we  present  in  detail  our  proposed  power  and  code¬ 
channel  allocation  solution.  The  performance  of  the  proposed 
scheme  is  evaluated  through  simulations  in  Section  IV.  A  few 
concluding  remarks  are  drawn  in  Section  V. 

II.  System  Model  and  Problem  Statement 

We  consider  a  primary  CDMA  system  with  processing  gain 
(code  sequence  length)  L,  K  primary  transmitters  P'I\,  i  = 
1,2, ...  ,K,  and  a  primary  receiver  PR  (for  example,  K 
uplink  transmissions  by  users  PTi,i  =  1.  2, ....  K ,  to  base 
station  PR).  We  also  consider  a  potential  concurrent  sec¬ 
ondary  code-division  link  in  the  spectrum  band  of  the  primary 
system  between  a  secondary  transmitter  ST  and  receiver  SR 
(Fig.  1).  All  signals,  primary  and  secondary,  are  supposed  to 
propagate  over  flat-fading  channels  and  experience  additive 

2The  given  SINR  constraints  on  the  primary  signals  may  be  drawn/derived 
from  preset  QoS  requirements  for  the  primary  system. 
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Fig.  1.  Primary/secondary  CDMA  system  model  of  K  primary  transmitters 
PTi,i  =  1,2,...,  K ,  a  primary  receiver  PR,  and  a  secondary  transmitter- 
receiver  pair  ST,  SR.  All  paths  hi, ... ,  h^,qi,  •  •  • ,  Qk->  Qs  exhibit 
independent  (quasi-static)  Rayleigh  fading. 

white  Gaussian  noise.  We  denote  by  hi,  qi,  i  =  1,2,...,  K, 
the  path  coefficients  from  PTt  to  PR  and  SR,  respectively. 
The  path  coefficients  from  ST  to  PR,  SR  are  denoted  by 
hs  and  qs,  respectively.  All  path  coefficients  are  modeled  as 
Rayleigh  distributed  random  variables  that  are  independent 
across  user  signals  and  remain  constant  during  several  symbol 
intervals  (quasi-static  fading). 

After  carrier  demodulation,  chip-matched  filtering  and  sam¬ 
pling  at  the  chip  rate  over  the  duration  of  a  symbol  (bit)  period 
of  L  chips,  the  received  signal  at  the  primary  receiver  PR  can 
be  represented  as 

K 

r  =  ^2  \[EihiSibi  +  \/E~shsssba  +  np,  (1) 

i—  1 

while  the  secondary  signal  received  by  SR  is 

K 

y  =  VPiQiSibi  +  \/l^qsssbs  +  ns  (2) 
2=1 

where  Ei  >  o,  bi  £  {±1},  and  Si  £  MA  ||sj||  =  1,  denote  bit 
energy,  information  bit,  and  normalized  signature  vector  of  pri¬ 
mary  user  i,  i  =  1,  2, . . . ,  K,  respectively;  Es  >  0,  bs  £  {±1}, 
and  ss  £  Rl,  ||ss||  =  1,  denote  the  bit  energy,  information  bit, 
and  normalized  signature  vector,  respectively,  of  the  secondary 
transmitter  ST;  np  and  ns  represent  additive  white  Gaussian 
noise  (AWGN)  at  PR  and  SR,  correspondingly,  independent 
from  each  other  with  0  mean  and  autocovariance  matrix  er2I. 

The  linear  filters  at  the  primary  and  secondary  receivers  that 
exhibit  maximum  output  SINR  [15]  can  be  found  to  be 

W maxSI N R,i  —  S^,  i  —  1,2,...,  K, 

WmaxSINR,s  =  C2R7lss, 


where  Rp  =  E{rrT},  Rs  =  E{ yyT},  ci,c2  >  0  (£{•} 
denotes  statistical  expectation  and  T  is  the  transpose  operator). 
The  output  SINR  at  PR  with  respect  to  the  signal  transmitted 
by  PTi,  SINRi,  is  given  in  (3)  followed  in  (4)  by  the  output 
SINR  at  SR,  SINRs,  where  Rp/,;  and  Rs/s  are  the  “exclude 
i”  and  “exclude  s”  input  data  autocorrelation  matrices  at  PR 
and  SR,  respectively,  defined  by 

K 

RpA  =  EkhlskSk  +  Esh2sssJ  +  er2I, 

k—l,k^i 

Rs/s  =  Ek1ksksl  +  cr2I. 

k= 1 

In  our  cognitive  radio  setup,  the  secondary  transmitter  has 
to  guarantee  the  SINR  QoS  of  all  primary  users.  In  this  spirit, 
our  objective  is  to  find  the  transmission  bit  energy  Es  and 
the  real-valued  normalized  signature  vector  ss  that  maximize 
SINRs  under  the  constraints  that  SINRi,  i  =  1,  2, . . . ,  K, 
are  all  above  a  certain  threshold  a  >  0,  i.e.  we  would  like  to 
identify  the  optimal  pair 

(Es,ss) opt  =  arg  max  Es s J R ss 

Ss>0,ssGlL 

subject  to  Eih2 sf  Rr);.  s i  >  a,  i  =  1, 2, . . . ,  K,  (5) 
sTs  —  1  E  <  K 

where  Emax  denotes  the  maximum  available/allowable  bit 
energy  for  the  secondary  user. 

The  optimization  task  of  maximizing  a  quadratic  objective 
function  (R“^  is  positive  definite)  subject  to  the  constraints  in 
(5)  is,  unfortunately,  a  non-convex  NP-hard  (in  L )  optimization 
problem  [16].  In  the  following  section,  we  delve  into  the  de¬ 
tails  of  the  problem  and  derive  a  novel  realizable  suboptimum 
solution. 


III.  Proposed  Cognitive  Secondary  Channel 
Design 

Using  the  matrix  inversion  lemma  [17]  on  R we  can 
express  the  key  quadratic  constraint  expression  sfRrhs i  in 
(5)  as 


^/i 


Si  = 


1  -  EihfsTHv  A 


*=1,2, 


,K,  (6) 


where  we  recall  that  Rp  =  /i{rry  }  is  the  autocorrelation 
matrix  of  the  whole  input  to  the  primary  receiver  PR.  Then, 
the  PR  SINR  constraints  in  (5)  become 


s'  R„  s,  > 


l-Eih2  +  aEih2 


A 

=  7 i, 


and  the  optimization  problem  can  be  rewritten  as 

{Es ,  ss)opt  =  arg  max  Es sf  RA  ss 

Bs>o,sseRi 

subject  to  sfRp  JSi  >  7j,  i  =  l,2,...,K, 
ss  ss  =  1,  Es  A  i7max. 


(7) 


(8) 


Using  the  matrix  inversion  lemma  on  Rp  1  this  time,  we  see 
that 


R,;1  =  R1 

P  p/s 


E‘h^p/eSs^tip/a 

1  +  EsPisjRrlsss 


(9) 
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SINRi  = 


Eil^axSINRAVElKbiSi)^} 


E{\WmaxSINR,i(^2k=l,kjH  V^fc/lfcSfc&fc  +  AE^hsssbs  +  np)|2} 
E{  |  ^ maxSINR,s  {y/E~sqsbs  ss)|2} 


=  E*hi*i  RP/,S 


SINR,  = 


E{\^TnaxSiNRJJ2k=i  VE-kqkskbk  +  ns)|2} 


=  Esq2s  sfR^Sg 


(3) 


(4) 


where  Rp/S  is  the  autocorrelation  matrix  of  the  input  to  the 
primary  receiver  PR  excluding  the  secondary  transmission, 

K  K 

Rp/S  =  E{(^2  \f~ElhiSibi  +  np)(^  +  np)T} 

i=  1  i=  1 

=  Z  £7</i*?s<8<’ + ct2l 
2—1 

Then,  inserting  (9)  in  (8)  we  can  express  the  optimization 
constraints  as  explicit  functions  of  the  code  sequence  of  the 
secondary  user  ss,  i.e. 


EshlsiR  }  ScsI’R  }  Si 
sTR-}  Si  >  s  p/s  s  p/s 


p'sl~  1  +  Esh2ss^R~jsss 


+  7i>  *  = 

(10) 


For  notational  simplicity,  define  the  L  x  L  matrix 


Bi  =  hi R-}ssiSjR-p}s  -  Piht R;/s  (11) 

where 

A  =  sjR-}sSi  -  7i,  i=  1, 2, ,  K.  (12) 

Then,  the  optimization  problem  in  (8)  can  be  rewritten  -for 
one  more  time-  as 


xopt  =  are  max  xTR  }  x 

subject  to  xTBiX  —  A  <  0, 

x  x  ^  Ernax 


i  =  l,2,...,K,  (13) 


where  x  is  the  amplitude-including  transmitted  signature 
vector  of  the  secondary  user,  x  =  y/E^ss.  We  notice  that 
for  (13)  to  be  solved  at  the  secondary  transmitter  ST,  the 
primary  receiver  PR  must  communicate  the  matrix  parameters 
B,  and  scalars  A>  *  =  1,2,...,  K .  Therefore,  no  explicit 
communication  of  the  primary  channel  codes  and  gains  is 
required  that  may  directly  compromise  the  privacy/security 
of  the  primary  system.  In  terms  of  the  computational  effort, 
however,  (i)  B,,  i  =  1,2,...,  A",  are  not  necessarily  pos¬ 
itive  semidefinite,  hence  the  problem  in  (13)  is  in  general  a 
non-convex  quadratically  constrained  quadratic  program  (non- 
convex  QCQP),  and  (ii)  the  complexity  of  a  solver  of  (13)  is 
exponential  in  the  dimension  L  (NP-hard  problem). 

To  circumvent  these  two  difficulties,  we  first  observe  that 
if  we  use  the  trace  property  of  matrices  U,V,  7>{UV}  = 
7>(VU),  we  are  able  to  represent  the  objective  function  in 
(13)  as 

xTr;/^  =  mr;/1sx}  (h) 

where  X  =  xxT.  Thus,  the  optimization  problem  in  (13)  takes 
the  new  equivalent  matrix  form 

Xopt  =  argmax  7>{R;AX} 
xeKLxl  s/s 

subject  to  7>{BiX}  <  A;  i=l,2,...,K,  (15) 
7>{X}  <  Emax,  X  >z  0,  ranfc(X)  =  1 


where  X  y  0  denotes  that  the  matrix  X  is  positive  semidefi¬ 
nite. 

So  far,  we  have  shown  that  the  original  secondary  cognitive 
link  design  problem  in  (5)  is  equivalent  to  the  one  in  (8),  (13), 
and  finally  (15),  and  is  non-convex  NP-hard.  To  effectively 
attack  the  problem  anyway,  we  now  propose  to  relax  the 
rank  constraint  in  (15)  and  proceed  by  solving  the  following 
problem  instead, 

X'  =  argmax  7>{RA  X} 

xeKIxI  s  8 

subject  to  7>-{BjX}  <  A >  i  =  1,  2, . . . ,  AT,  (16) 
MX}  <  Emax,  X  >  0. 

Then,  (16)  is  a  convex  polynomial-complexity  problem  that 
can  be  solved  using  semidefinite  programming.  Strictly  speak¬ 
ing,  we  can  solve  (16)  in  polynomial  time  within  an  error 
e  >  0  from  its  value  at  the  optimum  point  X'.  More  specifi¬ 
cally,  let  f0  =  MR;/sX}|x=x',  i.e.  f0  is  the  optimum  value 
of  the  constrained  (affine)  objective  function  in  (16).  Then, 
for  any  given  e  >  0,  semidefinite  programming  guarantees 
that  we  can  converge  in  polynomial  time  (polynomial  in  the 
input  size  L  and  in  the  error  requirement  function  log  1  /e) 
to  a  solution  that  lies  in  (fa  —  e,f0)  [18].  In  this  paper,  for 
the  semidefinite  programming  problem  in  (16),  we  propose 
to  use  a  primal-dual  interior-point  method  [19].  In  particular, 
we  consider  the  problem  in  (16)  as  the  primal  optimization 
problem,  we  create  a  differently  parameterized  equivalent  dual 
problem,  and  then  solve  both  problems  iteratively  in  a  coupled 
fashion.  Then,  each  iteration  can  be  implemented  in  0(L3) 
and  the  algorithm  converges  after  (){L  log  1  /e)  iterations  to 
the  matrix  X"  that  makes  the  objective  function  7r{R“^X| 
attain  a  value  within  (/Q  —  e,  f0).  The  proposed  method  is 
outlined  in  Fig.  2.  We  note  that  relaxing  the  rank  constraint 
of  the  non-convex  NP-hard  problem  in  (15)  we  created  the 
convex  optimization  problem  in  (16)  that  can  be  solved  in 
0(A4  log  1/e)  time  (by  semidefinite  programming  methods  as 
described  in  Fig.  2).  Of  course,  because  of  the  constraint  re¬ 
laxation  itself  the  objective  function  evaluated  at  the  optimum 
point  X'  in  (16)  is  just  an  upper  bound  on  the  value  of  the 
objective  function  evaluated  at  the  optimum  point  of  interest 
Xopt  of  (15),  MR7/6Xopt}  <  Ji-{R;/sX'}.  Moreover,  X' 
is  not  available  exactly  either  and  instead  we  have  X"  with 

MR^x"}  e  (MR/Ax'}  - e, MR;/Sx'}). 

To  summarize  our  developments  so  far  for  the  cognitive 
design  of  a  code-division  secondary  link,  first,  for  the  given 
primary  SINR-QoS  threshold  a  >  0,  we  test  whether  A. 
i  =  1,2,... AT,  in  (12)  are  all  greater  than  zero.  If  this 
is  not  true,  then  the  SINR-QoS  constraints  for  the  primary 
users  cannot  be  met  and  outright  no  secondary  transmission 
is  allowed  (see  flow-chart  in  Fig.  3).  Otherwise,  we  run  the 
procedure  of  Fig.  2  which  returns  matrix  X".  If  the  rank  of  X" 
is  1  with  eigenvalue,  eigenvector  pair  Ai,ai,  then  we  already 
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Formulate  the  pair  of  primal  and  dual  semidefinite  programming  problems: 

Primal  Dual 

max  Tr{R_)X}  min  brv 

XgRi-Ki-  1  Vs  J  velRK+i.ZeR1-*^ 

subject  to  7>{B^X}  <  /%,  i  =  1,  2, . . . ,  K,  subject  to  +  wjf+il  =  Ry*  +  Z 

Tr{X}  <  Emax,  X  ^  o  z  t  o,  Vi  >  0,  i  =  1, 2, . . . ,  K  +  1 

where  b  =  \J3 . . . ,  Pk;  Emax]  and  Vj  is  the  ith  element  of  v. 

(i)  Initialize  (X,  v,  Z)  <-  (1, 0, 1). 

(ii)  Calculate  AX,  Av,  AZ  by  method  in  Appendix. 

(iii)  Replace  AX  <-  i(AX  +  AXT). 

(iv)  Choose  0  <  r  <  1  and  select  steplengths  ni,  n 2 

nj  =  min(l,T  •  ni),  ni  =  sup{f?i  £l  :X  +  fj.AX  b  0}. 
ri2  =  min(l,  r  •  T12),  »2  =  sup{^2  €  M  :  Z  +  I2A Z  >z  0}. 

(v)  Update 

X  <— X  +  niAX, 
v  <—  v  +  n2Av, 

Z  <—  Z  +  n2AZ. 

(vi)  Repeat  (ii)-(v)  until  X  feasible  solution  for  primal  problem,  (v,  Z)  feasible  solution  for  dual  problem, 
and  the  difference  of  primal  and  dual  objective  values  lies  within  (— e,  e). 

(vii)  Return  X"  <—  X. 


Fig.  2.  Proposed  interior-point  algorithm  for  solving  (16). 

have  our  secondary  link  design  with  signature  ss  =  ai  and 
transmission  amplitude  Es  =  Ai.  If  the  rank  of  X"  is  not  1, 
further  work  is  needed  as  described  below. 

When  X'  of  (16)  (or  in  practice  X"  returned  by  Fig.  2) 
happens  to  be  of  rank  1  with  eigenvalue,  eigenvector  pair 
Ai,ai,  then  X'  =  Xopt  in  (15)  and  x°pt  =  vTiai  in  (13). 
Otherwise,  there  is  no  direct  path  from  X'  of  (16)  to  xopt 
in  (13).  In  this  case,  we  may  simply  consider  changing  the 
search  for  an  optimal  vector  in  (13)  to  a  search  for  an  optimal 
probability  density  function  (pdf)  of  vectors  that  maximizes 
the  average  objective  function  subject  to  average  constraints, 
i.e. 

/°pt(x)  =  argmax  £{xrR;/sx} 
fm  ' 

subject  to  i?{xTB,x}  <  j3 i,  i  =  1,  2, . . . ,  K,  (17) 

E{yi  x}  ^  Emax 

where  /(x)  denotes  the  probability  density  function  of  x.  This 
switch  to  a  statistical  optimization  problem  has  been  known 
as  the  “randomized  method”  in  semidefinite  programming 
literature  [18].  Using  the  commutative  property  between  trace 
and  expectation  operators,  the  pdf  optimization  problem  in 
(17)  takes  the  equivalent  form 

/°pt(x)  =  argmax  Tr{R^s  £{xxT}} 
fm  ' 

subject  to  7>{Bii?{xxT}}  <  /3j,  i  =  1, 2, . . . ,  K,  (18) 
7>{£,{xxt}}  <  Emax. 

We  can  show  that  fopt  (x)  is  in  fact  Gaussian  with  0  mean  and 
covariance  matrix  X',  fopt(x)  =  A/(0,X').  With  X"  from 
Fig.  2  as  a  close  approximation  of  X',  we  can  draw  now  a 
sequence  of  samples  x1;  x2, . . . ,  xp  from  Af(0,  X").  We  test 
all  of  them  for  “feasibility”  on  the  constraints  of  (13)  whether 
xjB iXp  <  fa,  Vi  =  1,2, ...  ,K,  and  xjxp  <  Emax,  p  = 
1,2 , ...  ,P,  and  among  the  feasible  vectors  (if  any)  we  choose 


Fig.  3.  Flow-chart  of  proposed  power  and  code  allocation  algorithm  for 
secondary  link. 

the  one,  say  x^°\  with  maximum  x3  Rt/s*-  °bjective  function 
value  (see  flow-chart  in  Fig.  3).  We  could  have  suggested  at 
this  time  a  cognitive  secondary  link  design  with  \/E^ss  = 
x*T .  Instead,  we  will  use  x!0J  as  an  initialization  point  to 
an  iterative  procedure  that  will  lead  to  a  much  improved  link 
design  vector.  The  iterative  procedure  is  developed  below  and 
its  performance  is  evaluated  by  simulation  studies  in  the  next 
section. 

First  we  express  Rs/s  as 

Rs/s  =  SS  ST  +  crT  (19) 

where  S  =  [si,  S2, . . . ,  Sk]  denotes  the  matrix  with 
columns  the  signatures  of  the  primary  users,  and  E  = 
diag(£^iq^,  •  •  ■ ,  Using  the  matrix  inversion 

lemma, 

r:/,  =  ^i-34S(e-'  +  Tsar's1-.  (20) 

Substitution  of  (20)  in  the  objective  function  of  (13)  leads  to 
xTR“/1  x  =  -^-xTx  -  ^txtQ  x  (21) 

s/s  cH  ct4 

where  Q  =  S(E-1  +  S)-1S3  .  In  (21),  the  first  term 

^■xTx  is  a  convex  function  while  the  second  term  -^xTQx 
is  a  concave  function  (which  implies  that  ^xTQ  x  is  convex). 
Based  on  the  first-order  conditions  of  convex  functions  [20], 
we  have 

xTx  >  2x(0)tx  -  x(0)Tx(0)  (22) 

where  x!l))  denotes  an  initial  feasible  vector.  Then,  we  com¬ 
bine  (21)  and  (22)  and  form  an  optimization  problem  that 
maximizes  the  following  concave  function 

Ax(0)Tx  _  3_xtq  x  _  _Ix(0)Tx(0)  (23) 

crz  a4  az 
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that  leads  to  a  suboptimum  solution  for  our  original  problem  in 
(13).  To  maximize  (23)  in  view  of  our  constraints  in  (13),  we 
restrict  all  non-convex  constraints  into  convex  sets  (lineariza¬ 
tion).  In  particular,  we  consider  the  non-convex  constraints 

xtB,x  —  fk<  0,  i  €  lnc ,  (24) 

where  Inc  denotes  the  set  of  all  indices  i  G  {1,2,...,  A"}  for 
which  xtB,x  is  a  non-convex  function.  Then,  we  decompose 
the  matrix  B,  into  its  positive  and  negative  parts 

B,  =  B+  -  B-  (25) 

where  B+  =  hi s?; sf It”/*  and  B~  =  Pih2aR~f3  are  all 
positive  semidefinite.  Therefore,  the  original  constraints  (24) 
can  be  written  as 


xtB+  x  -  fa  <  xTBj  x,  i  G  Xnc,  (26) 

where  both  sides  of  the  inequality  are  convex  quadratic 
functions.  Linearization  of  the  right-hand  side  of  (26)  around 
the  vector  x^0^  leads  to 


xtB+  x-ft  <  x(0)tB-  x(0)+2x(0)tB-  (x-x(°>),  i  G  lnc. 

(27) 

In  (27),  the  right-hand  side  is  an  affine  lower  bound  on  the 
original  function  xTB“  x.  It  is  thus  implied  that  the  resulting 
constraints  are  convex  and  more  conservative  than  the  original 
ones,  hence  the  feasible  set  of  the  linearized  problem  is  a 
convex  subset  of  the  original  feasible  set.  Thus,  by  linearizing 
the  concave  parts  of  all  constraints,  we  obtain  a  set  of  convex 
constraints  that  are  tighter  than  the  original  non-convex  ones. 
Now,  the  original  optimization  problem  takes  the  form 

x«  =  arg  max  4*(0)T*  -  4*TQx  - 

xel1  o1  er4  a1 

subject  to  xtB+  x  —  x^0,i  B~  (2x  —  x(0^)  —  Pi  <  0,  i  G  Xnc , 

xtB,x  -  pi<  0,  i  G  Inc ,  (28) 


xTx<  E 


max 


where  Xnc  =  {1,2,...,  A"}  —  Xnc ■  The  problem  in  (28) 
is  a  convex  QCQP  problem  and  can  be  solved  efficiently 
by  standard  convex  system  solvers  [21]  to  produce  a  new 
feasible  vector  x^1  L  The  objective  function  xTR“^x  in  (13) 
evaluated  at  x!|  j  takes  a  value  that  is  larger  than  or  equal  to  its 
value  at  x(°) .  Repeating  iteratively  the  linearization  procedure, 
we  can  obtain  a  sequence  of  feasible  vectors  x^°\  x*7\ 
x^2), . . . ,  x<77  with  non-decreasing  values  of  the  objective 
function  in  (13).  This  procedure  converges  after  very  few 
(eight  or  nine)  iterations  as  demonstrated  experimentally  in 
the  following  section. 


IV.  Simulation  Studies 

We  consider  a  primary  CDMA  system  with  signature  length 
(system  processing  gain)  L  =  16  and  K  synchronous  users. 
We  are  interested  in  establishing  a  secondary  code-division 
transmitter/receiver  pair  when  the  primary  system  is  fully 
loaded  to  overloaded,  say  AT  varies  from  16  to  20.  All 
signatures  for  primary  users  are  generated  from  a  minimum 
total-squared-correlation  optimal  binary  signature  set  which 
achieves  the  Karystinos-Pados  (KP)  bound  for  each  (A',  L) 
pair  of  values3  [22]-[24].  The  transmission  SNRs  of  the  AT 

3For  L  =  16,  when  K  <  L  the  KP-optimal  sequences  coincide  with  the 
familial'  Walsh-Hadamard  signature  codes. 


Fig.  4.  Secondary  transmission  percentage  as  a  function  of  the  number  of 
primary  users  under  Cases  rank{X.")  =  1  and  >  1  (the  study  includes  also 
the  code  assignment  scheme  in  [10]). 


primary  users  are  all  set  equal  to  =  15 dB,  i  =  1,  2, ...  AT; 
the  maximum  allowable  transmission  SNR  for  the  secondary 
link  is  set  to  E™gx  =  12 dB.  The  channel  coefficients  hi 
and  Qi,  i  =  1,2 (see  Fig.  1)  are  taken  to  be  the 
magnitude  of  independant  complex  Gaussian  random  variables 
with  mean  0  and  variance  4;  the  same  holds  true  for  hs  and 
qs.  The  receiver  SINR  threshold  for  primary  users  is  set  to 
a  =  2 dB  which  corresponds  to  an  average  raw  bit-error-rate 
(BER)  at  the  output  of  the  maximum  SINR  linear  filter  receiver 
of  about  10-1.  Ten  thousand  (10,000)  system/secondary-line 
optimization  experiments  are  run  under  the  described  (quasi¬ 
static)  flat  fading  conditions.  When  random  vector  drawing  is 
necessitated  by  the  flow-chart  in  Fig.  3,  P  =  50  test  vector 
points  are  generated. 

In  Fig.  4,  we  plot  as  a  function  of  the  number  of  primary 
system  users  AT  the  percentage  of  time  that  secondary  trans¬ 
mission  is  enabled  directly  under  the  case  rank(X.”)  =  1  or  by 
the  iterative  linearized  optimizer  as  well  as  the  “Interference- 
Minimizing-Code- Assignment”  (IMCA)  scheme  in  [10].  We 
observe  that  significant  opportunity  exists  for  cognitive  sec¬ 
ondary  transmission  when  the  primary  system  is  fully  loaded 
(K  =  L  =  16).  As  we  expect,  Fig.  4  shows  that  the  frequency 
of  secondary  transmissions  reduces  as  the  primary  system  load 
increases.  We  observe  also  that  our  proposed  scheme  offers 
more  opportunities  for  cognitive  secondary  transmission  than 
[10]. 

In  Fig.  5,  we  test  the  quality  of  the  secondary  transmission 
line  (the  pre-set  SINR-QoS  of  the  primary  system  is  -of 
course-  guaranteed  by  the  algorithmic  procedure)  and  the 
significance  of  the  iterative  linearized  optimizer  in  (28)  (see 
flow-chart  in  Fig.  3).  We  fix  the  primary  system  load  K  =  16 
(fully  loaded)  and  plot  the  secondary  receiver  average  SINR 
for  the  experimental  instants  of  ranfc(X")  >  1  as  a  function 
of  the  iteration  of  the  optimizer  initialized  at  the  point/design 
x(°>  that  is  the  best  out  of  P  =  50  samples  drawn  from  the 
A/”(0,X")  pdf.  It  is  pleasing  to  observe  that  eight  or  nine 
iterations  are  enough  for  effective  convergence. 
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Iteration  Number 


Fig.  5.  Secondary  receiver  SINR  as  a  function  of  the  iteration  step  of 
the  linearized  optimizer  in  (28)  initialized  at  the  best  feasible  sample  out  of 
P  =  50  drawings  from  the  J\f(0,  X")  pdf. 


Channel  Realization  Sequence 

Fig.  6.  Instantaneous  output  SINR  of  a  primary  signal  against  SINR-QoS 
threshold  a  (thick  line)  and  instantaneous  output  SINR  of  secondary  signal. 


Finally,  in  Fig.  6,  to  gain  visual  insight  into  the  operation 
of  the  primary/secondary  system  we  plot  the  instantaneous 
receiver  SINR  of  a  primary  signal  and  the  secondary  signal 
for  the  case  K  =  17  over  an  experimental  data  record  se¬ 
quence  of  1000  Rayleigh  fading  channel  realizations.  Missing 
secondary  signal  SINR  values  indicate  the  instances  when  no 
secondary  transmission  was  allowed.  The  proposed  scheme 
almost  doubled  the  occurrences  of  secondary  transmission 
compared  to  (10].  When  secondary  transmissions  do  occur 
for  both  schemes,  the  joint  power  and  sequence  optimization 
executed  by  the  proposed  scheme  results  in  superior  SINR 
performance  for  the  secondary  receiver  over  [10]. 

V.  Conclusions 

We  considered  the  general  problem  of  establishing  a  sec¬ 
ondary  code-division  line  alongside  a  primary  code-division 
multiple-access  system.  We  formulated  the  problem  as  the 
search  for  the  secondary  amplitude,  code  transmission  pair 
that  maximizes  the  secondary  line  output  SINR  subject  to 
the  condition  that  all  primary  signal  output  SINR  values  are 
maintained  above  a  given  SINR-QoS  threshold  value.  Regret¬ 
fully,  in  a  common  Rayleigh  fading  wireless  environment,  the 
formulated  constrained  optimization  problem  is  non-convex 
and  NP-hard  in  the  code  vector  dimension. 

Nevertheless,  in  pursuit  of  a  computationally  manageable 
and  performance-wise  appealing  suboptimal  solution,  we  first 
converted  the  amplitude/code-vector  optimization  problem  to 
an  equivalent  matrix  optimization  problem  under  a  rank- 
1  constraint.  Disregarding  (“relaxing”  in  formal  language) 
the  rank-1  constraint  makes  the  problem  amenable  to  an 
“easy”  polynomial-cost  semidefinite  programming  solution. 
When  luckily,  a  rank-1  matrix  happens  to  be  returned,  optimal 
secondary-line  design  is  achieved.  For  the  common  case  of 
a  higher  rank,  an  iterative  linearized  polynomial-cost  convex 
optimizer  is  developed  with  much  appealing  (yet  suboptimal) 
amplitude/code-vector  design  solutions  after  a  few  iterations. 

Cognitive  code-division  radios  combine  in  principle  the 
bandwidth  efficiency  characteristics  of  cognitive  operation  and 


code-division  multiple  accessing  and  are  expected  to  find 
a  place  in  future  communication  systems.  To  that  extend, 
the  developments  presented  in  this  paper  constitute  an  early 
contribution  that  can  be  helpful  in  benchmarking  future  efforts. 


Appendix  A 

Calculation  of  AX,  Av,  AZ  of  Fig.  2 

Let  vec{-}  denote  column-by-column  matrix  vectorization 
and  mat{-}  the  exact  inverse  operation.  Choose  0  <  <5  <  1 
and  define  /r  =  (>  r'{xz}  _ 

In  Fig.  2,  AX,  Av,  AZ  are  obtained  by  solving  the 
following  linear  system 


0  Bt  I 

vec{AX} 

uec{Tx} 

BOO 

Av 

= 

f  2 

EOF 

vec{  AZ} 

_  uec{T3} 

where 


(vec{Hi})T 


(29) 


&{K+l)xL2  — 


(vec{  BA'})T 
(vec{I})T 


tilxl  =  Rs/1s+Z-TOaf{B'rv},t2(K+1)xl  =  b— Bvec{X}, 

T3lxL  =/zI-XZ,El2xL2  =Z<g>I,  andFi2xi2  =I(g)X 
((££)  denotes  the  standard  Kronecker  product).  Applying  Gauss 
elimination,  the  solution  is 


Av  =  (BE_1FB)_1(t2  +  BE_1(F«ec{Ti}  -  ttec{T3})), 

(30) 

i>ec{AX}  =  — E_1(F(uec{T1}  -  BTAv)  -  t;ec{T3}), 

(31) 

vec{AZj  =  vec{T{\  —  BTAv.  (32) 
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